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| The two stars that form the subject of this paper are both short-period double-lined 

eclipsing binaries. They have non-circular orbits despite their short periods, both of which 

are very close to integral numbers of sidereal days (10.0015 and 9.0127 sidereal days 

respectively), making it difficult to achieve uniform phase coverage of the orbits from a 

■ single site in the short term. 

Although the HD type of HD 152028 is GO, the integrated spectral type of the system 

j^. ■ must actually be much earlier: the (B — V) colour index is only about m .37 and the Hip- 
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parcos parallax indicates an integrated absolute magnitude as bright as +l m .4. A published 
photometric investigation suggests that the primary star exhibits 8 Scuti pulsations, with 
a period of 0.1138 days. That period is not present in the radial velocities, but we have 
identified a comparable periodicity in the initially excessive residuals (a ~ 2 km s -1 ) in the 

Oh 1 

radial velocities of the primary star: there is an asymmetrical pulsational velocity curve with 
a semi-amplitude of ~3 km s and a period of 0.1178 days. The disagreement between 
that period and the published photometric one is probably because the latter is mistaken. 

It is argued from the luminosity of the system, the eclipse light-curves in V and 
(B — V), the lengths of the eclipse chords, and the rotational velocity of the primary 
star, that the eclipses are practically central (a total eclipse and a transit) and that the 
component stars have the same colour but luminosities that differ by l m .0. 

HDE 284195 was not observed by Hipparcos , but its HDE type of GO is in reasonable 
agreement with its colour and the nature of its radial- velocity traces. The rotations of both 
stars appear to be pseudo-synchronized to the orbit. 

Since this is 'only' a radial- velocity study and not a comprehensive discussion of the 
objects, the orbital inclinations are not formally determined. They must be very high, 
however (otherwise the systems would not eclipse), and there is specific evidence that the 
inclination of HD 152028 is very close to 90°, so it is permissible to assume that the factors 
sin^ i in the masses are close to unity. In that case the masses are scarcely above the 
minimum values, which in the case of HD 152028 are 1.78 and 1.42 Mq and in that of 
HDE 284195 are 1.10 and 1.01 M , with uncertainties < 1%. 
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Introduction 

This paper presents orbits for two ninth-magnitude binary stars to which attention has been called by the 
discovery of eclipses. Both of them have proved to be spectroscopically double-lined, so their orbital elements 
possess particular significance by giving the masses directly. Another characteristic, rather troublesome, that 
the systems share is the close approximation of the period of each to an integral number of sidereal days. 

HD 152028 is in Draco, about 1 /n of the way from ( Draconis towards /? Ursae Minoris; it was classified 
as type GO in the Henry Draper Catalogue. It had attracted no attention whatsoever from astronomers 
until Hipparcos 1 discovered its eclipses, to which it attributed an incorrect period of 16.960 days; the normal 
magnitude was found to be 8 m .76. The Hipparcos discovery resulted in the assignment 2 of the variable-star 
designation GK Draconis. Very recently Dallaporta et al? have published a photometric investigation that 
has identified the true period as 9.9742 days and has demonstrated that there are primary and secondary 
eclipses of virtually equal depths, m .36, but somewhat different widths, 0.050 and 0.039 of the orbital period. 
In addition there is 'photometric noise' that becomes less conspicuous during the secondary minimum and 
is accordingly attributed to the star that is being eclipsed at that phase; a 2.7-hour periodicity is identified, 
and the suggestion is made that the relevant component of the binary is a 5 Scuti variable. The Hipparcos 
parallax is 0".00337 ± 0".00069, corresponding to a true distance modulus of about 7 m .36 ± m .5; even 
without allowance for interstellar absorption, the absolute magnitude must be as bright as +l m .4 with the 
same uncertainty of m .5 — much brighter than would be expected for a system with an integrated spectral 
type of GO if it were a normal main-sequence pair. Hipparcos gives the (B — V) colour index as m .376, while 
Dallaporta et al? give it as m .35; it is clearly too blue for a GO type. 

HDE 284195 (BD +21° 605) is in Taurus, roughly midway between the Pleiades and the centre of the 
Hyades; at ninth magnitude, it was marginally too faint for inclusion in the original Henry Draper Catalogue 
but was picked up in the second volume of the Henry Draper Extension 4 ', where it was classed GO. Like 
HD 152028, it came to attention only when it was discovered to be an eclipsing system; the discovery was 
made by Kaiser 5 , who initially gave its period as 3.176 days. Soon afterwards, Kaiser and five collaborators 
issued a correction to the period, explaining 6 that the orbit must be non-circular, since the secondary eclipse 
did not come at phase .5: the originally-supposed period of 3.176 days was actually the interval between 
the secondary eclipse and the ensuing primary one, the new value for the period being 4.49407 days. Later 
still, Kaiser & Frey 7 , acting on information provided privately by L. Marshall in 1996, further revised the 
period by doubling it, giving a true value of 8.988476 ± 0.000007 days. Meanwhile, Kazarovets & Samus 8 
had gazetted HDE 284195 as a variable star with the designation V1094 Tauri. 

Despite all the interest in its photometry, it seems that no actual magnitudes of HDE 284195 have 
been published. The star was not on the Hipparcos programme, but was observed as a Tycho star; the 
new reductions in the Tycho 2 catalogue 9 give V T = 9 m .052 ± m .024, B T = 9 m .833 ± m .032. According to 
information to be found in the Hipparcos introductory volume (ref. 1, 1, p. 57) those magnitudes transform 
to V = 8 m .98, (B — V) = m .66, evidently with uncertainties of a few hundredths; a recalibration by Besscll 10 
of the Hipparcos and Tycho photometric system would suggest values of 8 m .97 and m .70. 

Radial velocities and orbits 

No radial velocities are known to have been published for either star, although Dallaporta et al. 3 men- 
tioned that spectroscopic observations of HD 152028 were in progress and the advice that Kaiser & Frey 7 
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ascribed to Marshall concerning HDE 284195 stemmed from radial-velocity observations that seem not to 
have been published. 

Both systems have recently been under observation with the Coravel radial-velocity spectrometer at 
the coude focus of the Cambridge 36-inch reflector. Observations of HD 152028 were begun in the 2001 
observing season and promptly demonstrated that the 16.960-day period announced by Hipparcos was in 
error. That is, as a matter of fact, no longer particularly noteworthy, as it has proved to be the case for 
quite a number of stars noted by Hipparcos as eclipsing binaries; this is the third example in this Magazine, 
after HD 44192 (V454 Aur) 11 and HD 98116 (FO Leo) 12 , while others, including HD 175900 (HP Dra) 13 
and HD 183361 (V2080 Cyg) 13 have been documented elsewhere. The period-finding routine used by the 
Hipparcos authors, who were of course obliged to deal with their photometric material on a wholesale basis, 
seems to have worked satisfactorily in cases where the variations were quasi-continuous, but did not cope well 
with eclipsing systems. In such cases, most of the photometric data are clustered at the level of maximum 
brightness, with relatively few outliers representing eclipses, and the opportunities for aliasing — particularly 
in the context of the very inappropriate typical Hipparcos 'window function' — were evidently a great deal 
higher than those authors appreciated. It has been pointed out previously (rcf. 11, p. 325) that radial- velocity 
observations, where the quantity being observed varies continuously instead of only on rare excursions, offer 
a much more favourable basis for establishing periods. Accordingly it did not take long to recognize from the 
radial-velocity data that the true period of HD 152028 is very close to 10 days. It is only fair to remark here 
that, although radial velocities have the advantage over photometry in the initial discovery of a periodicity, 
the position is reversed once the approximate period is known: whereas radial velocities change gradually over 
the whole period, the suddenness of the onset and cessation of eclipses means that the whole photometric 
range is spanned in a small fraction of the period and gives photometry an advantage by a factor typically 
of 10-100 times in determining the exact value. 

The HD 152028 system has continued under radial-velocity observation up to the time of writing of this 
paper, the total number of measurements available for the determination of the orbit being 50. They are 
set out in Table I. Although the orbital period is extremely close to the integral number of 10 sidereal days 
(= 9.9727 solar days), so close that the object comes to the meridian at virtually the same set of phases 
year after year, it has been possible to improve the phase distribution of the radial velocities by observing at 
different hour angles. That is facilitated by the high declination (+68°) of the object; HD 152028 is therefore 
circumpolar as seen from Cambridge (and indeed from latitudes all the way down to the Tropic of Cancer), 
and so can in principle be observed at any hour angle. In actual practice, however, it can not be observed 
with the Cambridge system at extreme hour angles because the coude beam of the 36-inch reflector travels 
up the polar axis and the Coravel occupies the space to the north of the telescope, obstructing the view of 
the sky below the North Celestial Pole. 

In almost all observations the velocities of both components have been found, usually but not invariably 
from a single trace that includes both 'dips', but the dip given by the secondary is extremely weak and often 
challenges the limits of what one can attempt to measure with the Coravel, cf. Fig 1. Even the dip that 
represents the primary is weak. It is obvious not only from the character of the radial- velocity traces but also 
from the absolute magnitude and from the colour index of the system (remarked upon in the Introduction 
above) that the integrated spectral type must be a great deal earlier than the GO given in the Henry Draper 
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Catalogue] that is borne out by the finding of Dallaporta et al? that the primary is itself an intrinsic variable 
of the 5 Scuti type, implying a spectral class of late A or early F. 

This is the point at which it is convenient for the writers to explain the origin of their collaboration. 
A paper only somewhat analogous to this one was drafted by the first-named author, submitted, and refereed 
by the second author. In the form in which it was initially submitted for publication, the r.m.s. residual of 
the velocities of the primary star of HD 152028 was as large as 2.0 km s _1 — much worse than might have 
been expected by comparison, for example, with the sub-km s -1 residuals for even the secondary component 
of HDE 284195. The residuals were not related to the published 3 5 Scuti periodicity of 0.1137601 days. In the 
course of refereeing the paper, H. M. J. B. identified in the residuals a somewhat different periodicity of 8.49 
cycles per day (period 0.1178 days). That discovery led to a collaboration that has resulted in a paper that 
is much more substantial and informative than the original one. 

We have been able to apply to the radial-velocity data a program intended for the solution of orbits of 
triple stellar systems, to solve simultaneously for the elements of the binary orbit and for the pulsation of 
the primary star. In doing so, we have necessarily treated the pulsational velocity curve as if it could be 
represented by a Keplerian orbit, but within the uncertainties of its determination that has proved to be 
a fully acceptable approximation. Much of the excess 'noise' that at first existed in the radial velocities of 
the primary star has disappeared as a result of the identification of the 0.1178-day period; the remaining 
r.m.s. residual is only 0.70 km s _1 . Even so, there is evidence (the presentation of which is deferred to 
the Discussion section below) that the pulsation is not strictly repetitive, so it is likely that its complete 
characterization would allow a modest further reduction in the primary star's radial-velocity residuals. 

Here, the principal concern is with the primary-secondary orbit, which is plotted in Fig. 2. The velocities 
that are plotted there for the primary star are not the directly observed values but have been corrected for the 
respective computed pulsational contributions. Those contributions are explicitly listed in the seventh column 
of Table I; the previous column tabulates the pulsational phase. To obtain near-equality of the weighted vari- 
ances of the velocities of the two components, it has been necessary to attribute to the secondary's velocities 
a weight of only y 2 5 ■ Table II presents the orbital elements of the binary system, and the characteristics of 
the short-period oscillation of the primary star in terms of the equivalent Keplerian elements. 

TABLE II 

Orbital and pulsational elements for HD 152028 (GK Draconis) 



Element 


Orbit 


Pulsation 


P (days) 


9.97380 


± 


0.00019 


0.1177753 ± 0.0000005 


T (MJD) 


52558.36 


± 


0.05 


52558.225 ± 0.004 


7 (km s _1 ) 


+2.23 


± 


0.11 




K x (kms- 1 ) 


64.97 


± 


0.13 


2.62 ± 0.17 


K 2 (kms- 1 ) 


81.14 


± 


0.6 




9 


1.249 


± 


0.010 




e 


0.0815 


± 


0.0021 


0.26 ± 0.06 


uj (degrees) 


82.0 


± 


1.7 


96 ± 13 


aisini (Gm) 


8.881 


± 


0.018 


0.0041 ± 0.0003 


a2 sin % (Gm) 


11.09 


± 


0.09 






0.281 


± 


0.002 




f(m 2 ) (M ) 


0.548 


± 


0.013 




misin 3 ? (M Q ) 


1.776 


± 


0.021 




m 2 sm3 * (Ms) 


1.422 


± 


0.013 




R.m.s. residual (kms -1 ) 


0.70 (Primary) 






3.5 (Secondary) 





4 



Although slightly fainter and of nominally the same spectral type, HDE 284195 is an altogether easier 
object to observe (and to discuss!) than HD 152028. Fig. 3 gives an example of a radial-velocity trace, 
showing two dips that, while not of very generous depth, readily give velocities whose characteristic accuracy 
is better than 1 km s _1 . The Cambridge observations did not start until 2002 September, so all of the 44 
measurements that are listed in Table III are from the 2002/3 observing season. Again the orbital period 
proves to be very nearly an integral number of sidereal days (in this case 9 days) , but the coincidence is not 
as exact as in the case of HD 152028: the period is longer by about eighteen minutes, so in each successive 
calendar month any given phase is presented at an hour angle that is progressively displaced about an hour 
westwards. That has perversely meant that towards the end of the observing season, when observations were 
necessarily made at westerly hour angles, they fell at much the same set of phases as those made earlier in the 
season when the free choice of hour angle usually led to their being made near the meridian! Thus the orbit 
plot (Fig. 4) shows conspicuously bunched data points, although the determination of the orbit is scarcely 
impaired by the bunching. Eight of the measurements have not been useable in the solution of the orbit 
owing to their being taken when the object was single-lined and the velocities of the two components were 
too similar to one another to be separately determinable; in those cases the observations were reduced as if 
they were truly single and the result plotted with an open symbol in Fig. 4. The 36 remaining observations 
showed HDE 284195 as double-lined or at least (in the case of blends) could be reduced as such; the velocities 
given by the weaker secondary dip were less accurately measureable than those of the primary, and it was 
found necessary to attribute a weight of y 3 to them to bring the variances for the two components into 
approximate equality. On that basis the orbital elements are: 



p 


= 8.9881 ± 0.0005 days 


COl3 = 


MJD 52656.260 ± 0.007 


7 


= +4.59 ± 0.07 km s" 1 


di sin i = 


7.772 ± 0.015 Gm 


Ki 


= 65.30 ± 0.12 km s" 1 


a 2 sin i = 


8.448 ± 0.0.024 Gm 


K 2 


= 70.98 ± 0.20 km s" 1 


f{mi) 


0.2321 ± 0.0013 M Q 


q 


= 1.087 ± 0.004 (= mi/m 2 ) 


.f(m 2 ) 


0.2980 ± 0.0026 M Q 


e 


= 0.2697 ± 0.0018 


mi sin 3 i — 


1.099 ± 0.007 Mq 


UJ 


= 333.2 ± 0.3 degrees 


rri2 sin 3 i = 


1.011 ± 0.005 Mq 




R.m.s. residual = 0.46 km s" 


1 (primary), 0.83 km 


s -1 (secondary) 



The radial velocities of the two stars that constitute HD 284195 are computed to differ by nearly 
33 km s _1 at the conjunctions, so it would be quite possible (and interesting) to watch the signature of 
the eclipsing component grow at the expense of that of the eclipsed one at such times, as was in fact done in 
the case (described in Paper 160 11 ) of HD 44192. In the present case, however, the object was inaccessible to 
observations at the relevant phases; moreover, it has been concluded 14 , after a sober assessment of the scientific 
value of such observations, that they do little more than serve as an inefficient substitute for photometry, so 
no excuse is really warranted for their absence. Primary eclipses (those in which the secondary star obscures 
the primary, at phase .241) will be favourably presented for observers near the longitude (0°) of Cambridge 
throughout the next observing season, every ninth night from 2003 September 1.13 to 2004 March 16.86. 

In the case of HD 152028 the components' radial velocities differ by scarcely more than 2 km s _1 at 
the conjunctions, so there would be no particular interest in making radial-velocity observations at eclipses: 
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the dip would appear single-lined but the proportions of it attributable to the respective components would 
vary. One can determine retrospectively that the primary eclipse, which is centred at phase .526, must have 
been in progress at the time of the single-lined observation of 2001 June 25, and that egress from eclipse had 
probably not quite finished at the times of the observations of 2001 November 1, 2002 September 26, and 
perhaps 2002 October 6. 

Discussion — HD 152028 

In Fig. 5 we present the Hipparcos photometry folded on the true orbital period, superseding the diagram 
in ref. 1, vol. 12. Following photometric convention, the zero-point of the abscissae in the Figure is taken as the 
phase of one of the conjunctions in the spectroscopic orbit (the one at which the eclipse deemed by Dallaporta 
et al. 3 to be the primary one takes place). In order to get the phasing of the eclipses observed by Hipparcos to 
agree with those of the conjunctions in the recently observed orbit, it has been found necessary to refine the 
period from the 9.97380 ± 0.00019 days determined above to 9.9742 days; the latter value is exactly the one 
found photometrically by Dallaporta et al. (9.9742 ± 0.0001). The change, though small in absolute terms, 
is somewhat uncomfortable in relation to the formal standard deviation of the spectroscopically determined 
period. The eclipse light-curves here look more plausible than those in the Hipparcos reference, where one 
of them is a chain of points in practically a vertical line. 9.9742 solar days equal 10.0015 sidereal days, so 
from the point of view of an observer at a fixed longitude it takes 1000 days for the phasing of observing 
opportunities to change by 0.15 days; opportunities would not average out uniformly until a slippage of a 
complete day had occurred, implying a lapse of getting on for 7000 days or say 18 years. 

All the evidence points towards a spectral type much earlier than the HD class of GO for HD 152028. 
There is the colour index of 1 ' 3 about m .37, the trigonomctrically-derived 1 absolute magnitude of about 
+ l m .4, the minimum masses of 1.78 and 1.42 M Q found from the orbital elements above, and finally the 
assignment 3 of part of the photometric variations to 8 Scuti pulsations, which do not occur in G-type stars. 
On the basis of actual masses that are scarcely greater than the minimum ones, which must be the case 
since the fact that the system exhibits eclipses demonstrates that sin i ~ 1 , one could expect main-sequence 
types of about A9 and F3. Confusingly, however, Dallaporta et al? found that there is no perceptible change 
of {B — V) colour index during cither of the eclipses: the implication is that the components have identical 
colour indices despite their very unequal masses. 

It is the primary that Dallaporta et al. found to be the pulsating star, since the 'photometric noise' that 
they attributed to pulsations was considerably reduced when that star was in eclipse. The light-curves at 
both eclipses appear to have rather flat minima, as if they represented a total eclipse and a transit, but in 
both cases the relatively flat part at the bottom of the minimum is not truly horizontal but slopes down 
towards egress, so the flat-bottomed appearance may be a misleading effect of noise. Since the components 
seem to have identical colour indices they can be supposed also to have identical surface brightnesses. In that 
case the m .36 depths of both eclipses tell us that at each eclipse there is a loss of 28% of the total of the 
stars' visible surface areas. That would agree with the idea of the eclipses being total, with the stars having 
relative surface areas, and therefore luminosities, in the ratio 72 to 28, very close to one stellar magnitude. 
Their actual absolute V magnitudes would need to be +l m .8 and + 2 m .8 to constitute the +l m .4 found by 
Hipparcos for the system as a whole. 
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The relative transverse velocities of the component stars at the times of the eclipses are computed to be 
134 km s _1 at the conjunction where the primary is eclipsed (on the apastron side of the orbit — the eclipse 
called by Dallaporta et al. 3 "Min. I") and 157 km s _1 at the other. Dallaporta et al. say that those eclipses 
last for about 0.050 and 0.039 of the orbital period, i.e. 0.50 and 0.39 days, respectively. Simply by multiplying 
together the durations and velocities, we obtain the lengths of the eclipse chords: they are respectively 5.8 
and 5.3 Gm. It is physically impossible for the length of the eclipse chord on the apastron side of the orbit 
to be longer than that at the eclipse near periastron, since any departure of the orbital inclination from 90° 
must result in the apastron chord being shorter than the other. Dallaporta et ai.'s graphs of the eclipse 
photometry demonstrate that photometric noise seriously limits the accuracy with which the durations of 
the eclipses can be determined, so we probably do not need to worry about the apparent discrepancy, and for 
the purposes of the present discussion we adopt a mean eclipse chord of 5.55 Gm or 8.0 Rq for both eclipses. 
Even if the eclipses are central (orbital inclination 90° exactly), that means that the sum of the radii of the 
stars is 4.0 Rq, at first sight a surprisingly large value that goes a long way to reinforce the idea that the 
eclipses are practically central. 

On the basis (proposed above) that the stars have equal surface brightnesses and that the areas of their 
surfaces are in the ratio 72 to 28 or 2.57 to 1, their individual radii must be approximately 2.5 Rq for the 
primary and 1.5 Rq for the secondary. The latter radius is 10-20% larger than the tabular 15 value for a 
main-sequence star of the colour index of the HD 152028 system (F3V, m .37), but is actually near the 
lower bound of the distribution of radii with colour index in the graph shown by Andersen 16 of values well 
determined from eclipsing binaries, so we can view it with some equanimity. Neither the radius nor the colour 
of the primary star is appropriate to a main-sequence star of its mass, but we might best regard that star not 
as an anomalous object but as a somewhat evolved one, having begun its evolution towards the giant branch 
in the H-R diagram by moving towards the right from an initial position corresponding to a main-sequence 
type of about A9. Its present type might be estimated at F2 III-IV. 

We turn now to the evidence for pulsational instability of the primary star. Dallaporta et al. 3 were able 
to reduce the initially apparent 'photometric noise' by demonstrating in it a periodicity, which they attributed 
to S Scuti variability, of a fraction of a day. Similarly through the identification of a quite similar period 
in the radial velocities, a major reduction has been effected in what initially appeared to be 'radial-velocity 
noise'. 

There is an obvious embarrassment over the question of the pulsational period because, whereas the 
photomctrists have identified a period of 0.1136701 days (curiously enough stated to seven significant figures 
but attributed an uncertainty of as much as 0.0003 days), the period that suits the radial velocities is 
0.1177753 days with an uncertainty of only half a millionth of a day or some 40 milliseconds. We attribute 
the discrepancy to an error in the photometric period, which was established from measurements that were all 
taken within the same small range of orbital phases and therefore at epochs differing by multiples of 9.974 days. 
During one orbital period the number of pulsational periods according to the photometric value is 87.677, 
while according to the radial-velocity value it is 84.676. It may well be supposed that the photometric period 
could not be determined within one night to a precision sufficient to permit an unambiguous cycle count to be 
made to an epoch 10 days away (the actual intervals between the four nights utilized by Dallaporta et al. 3 in 
their pulsational analysis were 40, 50, and 10 days), and they made an error of three cycles per orbital period. 
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If they had not seemingly arbitrarily restricted themselves to photometry made on one particular night of 
the orbital cycle they would have picked up the mistake. We can demonstrate the merit of the radial- velocity 
period photometrically from the Hipparcos data: in Fig. 6 they are plotted on the photometric period and 
in Fig. 7 on the radial-velocity one. Only in Fig. 7 do the data show a systematic run with phase. It should 
be noted that the phasing of the photometric wave demonstrated in Fig. 7 is not secure: in the ~4000 days 
or 35,000 photometric periods since the Hipparcos epoch the 1-tr uncertainty in the period multiplies up to 
an interval of 0.017 days or 0.14 pulsational periods. 

The short-period radial- velocity variation was established from observations that were made on different 
nights and therefore in cycles far removed from one another. The observed variation must be slightly blurred 
because the radial-velocity integrations typically lasted nearly half an hour or 0.02 days, about one-sixth of 
the pulsational period. There seemed to be possible merit in observing the pulsation in 'real time', and by 
ignoring the secondary component of the binary system the integration times could be somewhat curtailed 
without loss of precision in the velocities of the primary. To that intent HD 152028 was observed continuously 
in a dedicated four-hour interval on the night of 2003 April 17/18. Only the region of velocity space that 
included the primary dip was scanned, in a series of fifteen 960-second integrations starting at 1017-second 
intervals, the latter being one-tenth of the period. The results are listed in Table IV. 



TABLE IV 

Additional radial-velocity observations of HD 152028 A 



Eel. Date (UT) 


Velocity 


Eel Date (UT) 


Velocity 


Eel. Date (UT) 


Velocity 




km s~ J 




km s^ 1 




km s~* 


2003 April 17.988 


+51.7 


2003 April 18.047 


+43.4 


2003 April 18.105 


+48.9 


18.000 


+52.3 


18.058 


+42.5 


18.117 


+49.0 


18.011 


+50.5 


18.070 


+45.1 


18.129 


+47.2 


18.023 


+47.7 


18.082 


+46.3 


18.141 


+44.2 


18.035 


+44.9 


18.094 


+47.4 


18.152 


+40.7 



In Figs. 8 and 9 the observed radial velocities are plotted against pulsational phase after the large 
contribution from the computed orbital velocity variation has been removed; the points corresponding to 
the observations listed in Table I and in Table IV are distinguished from one another by the use of different 
plotting symbols. In Fig. 8 the velocity curve is that derived from the observations in Table I alone and 
corresponds to the pulsational elements given in Table II. It is seen that the Table IV velocities do not fit 
it very well: they suggest a curve of considerably larger amplitude and possibly slightly shifted phase. The 
difference in amplitudes is much greater than could be explained by the somewhat sharper time-resolution 
of the new measurements. The curve in Fig. 9 is drawn to suit principally the Table IV points, although 
the others had to be retained in the solution at some level (they were given a weight of Vso), because they 
are necessary for the derivation of the orbital velocity curve which must be known before the pulsational 
contribution to the velocities can be isolated. The r.m.s. residual of the Table IV points from the curve 
drawn to fit them is only 0.41 km s _1 , but the curve definitely does not suit the other points as well as the 
one plotted in Fig. 8. (There is not a one-to-one correspondence between the positions of the points plotted 
in the two diagrams, owing to the re-computation of the orbital velocity variation with the greatly changed 
weighting system of Fig. 9.) 

This exercise tends to show that the pulsation responsible for the photometric and radial-velocity 
oscillations, while maintaining clockwork regularity in phase over an interval of years, does not maintain 
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its form very exactly. A tentative conclusion, suggested by analogy with other examples of pulsating stars, 
might be that there are additional periodicities with periods rather similar to the principal one, producing a 
beat phenomenon that is manifested in the apparent discrepancy between the velocity curves seen in Figs. 8 
and 9. It seems fairest to adopt — as has been done in Table II and is illustrated in Fig 8 — the curve that 
is assembled from observations made in random and well separated cycles, as being more representative of 
the general behaviour of the star than the curve determined at a single epoch from the observations listed in 
Table IV and plotted in Fig. 9. 

It is possible to use the relationship between period, luminosity, mass, and temperature (e.g. ref. 17, 
Equation 2) to obtain the bolometric magnitude of HD 152028 A as a function of the pulsation constant Q. 
We derive: 

My = 6.78 + 3.33 log Q, 

where we have assumed a bolometric correction of m . 1 1 . All the quantities that needed to be entered into the 
period-luminosity relationship to obtain the above result are very well constrained. If the period of 0.117775 
days that we have detected represents the fundamental radial mode, for which the corresponding value of Q 
is 17 0.033 days, we obtain My = +l m .85, in full agreement with the value obtained above from the Hipparcos 
photometry and the eclipse depth. 

The difficulty of measuring the radial velocity of the secondary star has seemed not to be the same on 
different occasions, as if the secondary too were of varying character, and certainly many of the residuals from 
the computed velocity curve are appalling and without precedent in this work. There is no clear periodicity 
in the residuals; although a period analysis seems to indicate a possible period close to 0.4 days we would not 
give much for its reality If the seeming variations and residuals stem simply from observational difficulties 
and errors resulting from the marginal observability of that star they are regretted, but the possibility is not 
absolutely ruled out. We note that the position of the star in the H R diagram is within the area populated 
by 7 Doradus variables, which have periods in the range 0.4-3 days. 

In the absence of adequate experience concerning stars such as those involved in HD 152028, no reliable 
deductions can be made from the observed strengths of the dips seen in radial- velocity traces. The widths of 
those dips may nevertheless be significant. The widths found from radial- velocity traces of the secondary star 
are so uncertain and scatter so widely that it would be unwise to try to make any deduction from them. The 
mean value for the primary, however, is quite accurately established; if interpreted purely as the effect of stellar 
rotation, it yields a vsini of 14 km s _1 . In view of the very high inclination that we know to characterize the 
orbit and that can reasonably be expected to apply also to the rotational axis, the same figure could be taken 
to represent the actual equatorial velocity. If we further assume the rotation to be pseudo-synchronized 18 
with the orbital revolution, the rotation period is 9.6 days and a stellar radius of 2.6 Rq is derived. It could 
be argued that the stellar pulsation may broaden the line profile and falsify the vsini value derived from it, 
so 2.6 Rq must represent a maximum value. The pulsational amplitude is, however, small compared with 
the apparent rotational velocity, so its contribution to the line-width is probably minimal. Moreover, the 
similarity of the derived radius to the one found from the length of the eclipse chord — a minimum value 
that seems difficult to refute — encourages a belief that the various considerations are converging towards a 
realistic model of the system, one that harmonizes, moreover, with the trigonometrical parallax. 
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Discussion — HDE 284195 

HDE 284195 seems to be a much more straightforward double-lined binary system than HD 152028. The 
disparity in equivalent widths of the two dips seen in radial- velocity traces (cf. Fig. 3) suggests a difference 
in V luminosity of m .5-0 m .6, corresponding to about three spectral sub-types; the mass difference of about 
9% agrees. The masses of 1.10 and 1.01 M given by the orbit above on the assumption (validated by the 
eclipses) that sin i ~ 1 would suggest the types to be F9 V and G2 V; the integrated colour index of m .66 (or 
according to the Bessell calibration 10 m .70) derived from the Tycho 2 photometry might urge us to adopt 
slightly later types, such as GO + G3 or even Gl + G4, but it is to be recalled that the standard error of 
that index is as much as m .04. 

As in the case of HD 152028, the observed projected rotational velocities vsini for the components of 
HDE 284195 may be considered as the actual equatorial velocities. The mean values are 9 and 7 km s _1 for 
the primary and secondary, respectively, with formal standard deviations that are smaller than the lower limit 
of 1 km s _1 that is ever claimed for rotational velocities determined from radial- velocity traces. If the radii 
of the stars are estimated at 1.1 and 1.0 Rq, respectively, the periods of their rotations are found to be 6.2 
and 7.2 days. At the orbital eccentricity of 0.27 found for HDE 284195, the pseudo-synchronous rotational 
period 18 is shorter than the orbital period by a factor of 1.49, making it 6.0 days, so we can be pretty certain 
that both stars are pseudo-synchronized. 

Note added in proof 

A preprint of a paper to be published in A&A by Zwitter et al. (a consortium whose membership overlaps that of Dallaporta 
et ai. 3 ), referring in part to HD 152028, has just appeared. It gives radial velocities that seem to be more scattered than those presented 
in the present paper and does not throw any light on the pulsational period. In fact a light-curve analysis is attempted on the basis of 
a plot (exactly like Fig. 5 here) of the raw Hipparcos and Tycho photometry, by way of demonstrating what might be obtainable from 
such data rather than in an effort to obtain the best possible values for the stellar parameters. 
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Fig. 1 

Radial-velocity trace of HD 152028 (GK Draconis), obtained with the Cambridge Coravel on 2003 March 1 and showing the 
unequal dips. 

Fig. 2 

The observed radial velocities of HD 152028 (GK Draconis) plotted as a function of phase, with the velocity curves corresponding 
to the adopted orbital elements drawn through them. Squares represent the velocities of the primary, circles those of the secondary 
The open diamonds indicate 'single-lined' measurements, where the velocities of the components were so similar to one another that 
they could not be assigned separately; they were not used in the solution of the orbit. The measurements of the primary star have been 
corrected for the pulsational radial- velocity variations on the basis of the averaged pulsational 'elements' listed in Table II. 

Fig. 3 

Radial- velocity trace of HDE 284195 (V1094 Tauri), obtained with the Cambridge Coravel on 2003 February 18. 

Fig. 4 

The observed radial velocities of HDE 284195 (V1094 Tauri) plotted as a function of phase, with the velocity curves corresponding 
to the adopted orbital elements drawn through them. Squares represent the velocities of the primary, circles those of the secondary. 
The open diamonds indicate 'single-lined' measurements, which were not taken into account in the solution of the orbit. 

Fig. 5 

The Hipparcos photometry of HD 152028, phased to the true orbital period of 9.9742 days. 

Fig. 6 

The Hipparcos photometry of HD 152028, phased to the pulsational period given by Dallaporta et a/A 

Fig. 7 

The Hipparcos photometry of HD 152028, phased to the pulsational period found in the radial velocities of the primary star. The 
fact that the photometry shows a clear variation with phase in this diagram (and not in Fig. 6) provides convincing evidence that this 
period is the correct one. The outlying low points (there are others that fall below the bottom of the diagram) are explained by the 
eclipses in the system. There would be obvious advantage in 'cleaning' the photometry of the pulsational variation documented by this 
Figure before using it in any discussion of the eclipse light-curves. 

Fig. 8 

The pulsational radial- velocity curve of the primary star in HD 152028, established from the observations given in Table I (plotted 
here as filled squares) by a simultaneous solution of the orbital and pulsational variations. The 'elements' of the latter, specified as if 
they were Keplcrian orbital elements, are included in Table II. The open star symbols represent radial- velocity observations made in a 
continuous run covering F/2 cycles of the pulsation on the night of 2003 April 17/18. They were not utilized at all in the solution of the 
orbital and pulsational elements; they are seen to confirm the variation in a general way, but to exhibit a distinctly greater amplitude. 

Fig. 9 

As Fig. 8, but this time the starred observations were given by far the greatest weight in the calculation of the orbital and 
pulsational elements. To assist the eye to appreciate which points the computed velocity curve is trying to match, the filled or open 
characters of the plotted points have been reversed in comparison with Fig. 8. The largcr-than-avcragc amplitude of the pulsation on 
the night when the 'star' points were observed is obvious. 
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TABLE I 



Radial-velocity observations 

Heliocentric Date HMJD Velocity 

Prim. Sec. 











km s 


1 

- 1 


km s 


1 

— 1 


2001 May 13 


.120 


52042 


.120 


-59, 


1 


+80, 


.7 


June 18 


.040 


078 


.040 


+67, 


,5 


-76. 


A 


23 


.034 


083 


.034 


-54, 


3 


+70. 


.9 


25 


.005 


085 


.005 




+2 


.0 




28 


.020 


088 


.020 


+65, 


,6 


-73 


.4 


July 15 


.994 


105 


,994 


+41. 


,0 


-43. 


.0 


25 


.971 


115 


.971 


+43, 


,8 


-53. 


.6 


26 


.962 


116 


.962 


+63, 


,2 


-68. 


,4 


28 


.030 


118 


.030 


+64, 


,0 


-77. 


,2 


30 


.034 


120 


.034 


-12, 


,5 


+13. 


.8 


Aug. 1. 


.042 


122 


.042 


-60, 


,7 


+78. 


,4 


1. 


,918 


122 


,918 


-49, 


,4 


+68, 


.9 


Nov. 1. 


.769 


214 


.769 


+ 13, 


,2 


-13 


,2 


2002 Apr. 27 


.094 


52391 


.094 


-60, 


,4 


+75. 


.1 


May 28 


.070 


422 


,070 


-51, 


,6 


+70. 


.8 


July 26 


.943 


481 


.943 


-54, 


,8 


+73. 


.0 


Sept. 2. 


.900 


519 


.900 


-47, 


,2 


+68. 


,0 


3 


.949 


520 


.949 


-61. 


,1 


+79. 


,1 


4. 


.965 


521 


.965 


-52, 


.6 


+68. 


.3 


12 


.833 


529 


.833 


-45, 


,9 


+63. 


.1 


26 


.892 


543 


.892 


+ 10, 


3 


-5, 


,0 


28 


.877 


545 


.877 


+61, 


,5 


-73, 


,5 


Oct. 6. 


.895 


553 


.895 


+9, 


5 


-7, 


,1 


12 


.869 


559 


.869 


-48, 


.8 


+68. 


,7 


23 


.800 


570 


,800 


-59, 


,4 


+83. 


,5 


Nov. 4. 


.784 


582 


.784 


-23, 


,7 


+36. 


.0 


Dec. 4. 


,741 


612 


.741 


-20, 


,1 


+35. 


,2 


9. 


.720 


617 


.720 


+30, 


,7 


-34, 


.0 


2003 Jan. 11 


.278 


52650 


.278 


-59, 


,6 


+86. 


.3 


16 


.250 


655 


,250 


+63, 


,0 


-69. 


,1 


28 


.237 


667 


,237 


+49, 


,6 


-53, 


,4 


Feb. 15 


.227 


685 


.227 


+57, 





-73, 


,4 


19 


.089 


689 


.089 


-38, 


,7 


+55. 


,0 


20 


.078 


690 


.078 


-59, 


.8 






21 


.206 


691 


.206 


-57, 


,7 


+82, 


.6 


22 


.089 


692 


.089 


-35, 


.0 


+46. 


,8 


Mar. 1. 


,223 


699 


.223 


-44, 


,6 


+58. 


3 


3 


,196 


701 


.196 


-54, 


,7 


+76. 


.5 


15 


.208 


713 


.208 




+5 


.3 




16 


.160 


714 


.160 


+34, 


,5 


-35. 


.3 


17 


.102 


715 


.102 


+58, 


,8 


-69. 


,7 


19 


.140 


717 


.140 


+44, 


,2 


-54, 


,1 


20 


.123 


718 


.123 




+4 


.6 




23 


.142 


721 


.142 


-57, 


.6 


+73, 


,2 


24 


.000 


722 


,000 


-35, 


,7 






Apr. 1. 


.045 


730 


.045 


-59, 


,5 


+85, 


.9 


7. 


,153 


736 


.153 


+70. 


.0 


-80. 


.7 


8. 


.069 


737 


.069 


+47. 


.6 


-50. 


.8 


16 


.108 


745 


,108 


+63, 


.6 


-71, 


,8 


17 


.031 


746 


.031 


+67, 


,8 


-77. 


,2 



of HD 152028 (GK Draconis) 

Phase V r pu i se (O - C) 
Orbit Pulse Prim. Sec. 

km s^ 1 km s' 1 km s~ 






.240 


0, 


,884 


+2 


.2 


+0, 


.6 


-1, 


.7 


3 


.841 


305 


,871 


+2 


.4 


-0, 


.1 


+0, 


.0 


4 


.342 


348. 


,276 


-2 


.0 


-0 


,5 


+1. 


,2 




,540 


365. 


,012 


-0 


.7 












,842 


390 


,608 


+1. 


.2 


-0, 


.8 


+2. 


.9 


6 


.644 


543. 


,224 


-2 


.4 


+1. 


,3 


+4. 


.6 


7 


.645 


627. 


,936 


+1. 


.4 


+0. 


,2 


-5. 


.9 




,744 


636 


,349 


-1 


.3 


+0, 


,3 


+6. 


.8 




,851 


645. 


,415 


-0 


.7 


+0. 


.8 


-2. 


.4 


8 


.052 


662 


,429 


-0 


.5 


+0 


,8 


-7. 


.1 




,253 


679. 


,481 





.0 


+1. 


,2 


-4. 


.0 




,341 


686 


,920 


+1 


.7 


+0, 


.9 


-1, 


.1 


17 


.550 


1466 


,805 


+2 


.5 


-0, 


.8 


-4, 


.0 


35 


.229 


2963. 


,932 


+1. 


.5 


-0. 


,2 


-6. 


.9 


38 


.335 


3226 


,943 


+1. 


.2 


+0. 


,4 


-0. 


,7 


44 


.338 


3735 


.309 


-1 


,7 


-0, 


,5 


+2. 


.3 


48 


.144 


4057. 


.589 


+ 1 


.1 


-0. 


,4 


+3 


,2 




.249 


4066. 


,498 


+0 


.2 


+0 


,7 


-3 


,4 




,351 


4075. 


,126 


-2 


.6 


+0. 


.1 


+0. 


,8 


49 


.139 


4141. 


,928 


+1. 


.5 


-0, 


.8 


-0. 


,2 


50 


.549 


4261. 


,304 


-1 


.8 


+1. 


.1 


+3, 


.7 




.748 


4278, 


,156 


-2 


.7 


-0, 


.6 


+2, 


.4 


51 


.552 


4346, 


,237 


-2 


.3 


-0, 


.2 


+2, 


.9 


52 


.151 


4396 


.956 


+0 


.9 


+0. 


.1 


+1. 


,4 


53 


.247 


4489. 


,770 


+2, 


.5 


+0. 


,2 


+1. 


.0 


54 


.448 


4591 


,526 


+0 


.5 


-0. 


.9 


+1. 


.9 


57 


.452 


4845. 


,883 


+2, 


.2 


-0. 


,2 


+2. 


.5 




.951 


4888. 


,153 


-2 


.7 


-0. 


,8 


+3. 


.6 


61 


.216 


5164. 


,596 


+ 1 


.1 


0, 


.0 


+5, 


,4 




.714 


5206. 


,814 


+2, 


,6 


+1. 


.1 


0, 


.0 


62 


.916 


5308. 


,589 


+ 1 


.1 


+0, 


.9 


+1. 


.1 


64 


.720 


5461. 


,341 


-1 


,4 


-2. 


.0 


-3 


.0 


65 


.107 


5494. 


,129 


-2 


.6 


-0, 


,1 


+5, 


.1 




.206 


5502. 


.528 


+0 


.5 


-0. 


,4 








.319 


5512. 


,107 


-2 


.5 


+0, 


,7 


+7. 


.7 




.408 


5519. 


.606 


+1. 


.2 


-0 


,5 


-2. 


,8 


66 


.123 


5580. 


,174 


-2 


,7 


-0. 


,2 


+1. 


,4 




.321 


5596. 


,930 


+1. 


.5 


-0 


,5 


+2. 


.0 


67 


.525 


5698. 


.922 


+1. 


.6 












.621 


5707. 


,005 


-0 


.5 


+0 


.1 


+3. 


.3 




.715 


5715. 


,002 


-0 


.4 


-0 


.3 


-0. 


,4 




.919 


5732, 


,303 


-1 


.8 


-0, 


.5 


-1, 


.1 


68 


.018 


5740, 


,649 


+ 1 


.6 












.321 


5766, 


,286 


-1 


.9 


0, 


.0 


-1, 


.4 




.407 


5773, 


,572 


+0 


.9 


-0. 


.6 






69 


.213 


5841. 


,877 


+2 


.3 


-1, 


.2 


+5. 


.2 




.826 


5893. 


,739 


+2, 


,3 


+0. 


.9 


-2 


,2 




.917 


5901 


,515 


+0 


.4 


+0. 


,2 


+3. 


.0 


70 


.724 


5969. 


,775 


+2 


.5 


+0. 


.0 


-0. 


,5 




.816 


5977. 


,611 


+1. 


.3 


-0 


.9 


+2 


.1 



TABLE III 



Radial-velocity observations of HDE 284195 (V1094 Tau) 

Heliocentric Date HMJD Velocity Phase (O — C) 

Prim. Sec. Prim. Sec. 

km s' 1 km s^ 1 km s' 1 km s' 1 



Sept. 28.096 


52545.096 


-44.9 


+57.8 


0.632 


0.0 


-0.6 




30.149 


547.149 


+3 


.7 


.861 






Oct. 


4.136 


551.136 


+2. 


.9 


1.304 








18.107 


565.107 


+3 


.4 


2.858 








19.097 


566.097 


+64.8 


-61.1 


.969 


0.0 


-0.2 




24.085 


571.085 


-40.2 


+53.9 


3.524 


+0.2 


+0.5 




28.144 


575.144 


+68.0 


-64.4 


.975 


-0.2 


+0.1 


Nov. 


2.127 


580.127 


-40.1 


+54.9 


4.530 


+0.8 


+0.9 




7.070 


585.070 


+80.7 


-78.2 


5.080 


0.0 


0.0 




11.035 


589.035 


-40.5 


+53.1 


.521 


-0.4 


-0.1 


Dec. 


4.974 


612.974 


+41.9 


-34.5 


8.184 


-0.2 


+1.6 




9.883 


617.883 


-38.2 


+49.3 


.730 


-0.8 


-1.0 




11.083 


619.083 


+5. 


.1 


.864 







Jan. 


5. 


,000 


52644 


.000 


-45. 


.3 


+58.4 


11.636 


-0 


,4 


0, 


.0 




5, 


.978 


644 


.978 


-35. 


,6 


+47.8 


.745 


-0 


.6 


+0. 


.2 




7, 


.000 


646 


.000 




+4, 


1 


.858 












7. 


.853 


646 


.853 


+56. 


,4 


-51.2 


.953 


+0 


.1 


+0, 


.5 




10. 


,023 


649 


.023 


+37. 


.6 


-31.3 


12.195 


-0 


.1 


+0. 


.2 




10. 


.975 


649 


,975 




+4. 


.1 


.301 












11. 


,766 


650 


.766 


-21. 


.8 


+32.0 


.389 


-1 


.2 


+0. 


.1 




11. 


.890 


650 


.890 


-22. 


.9 


+34.2 


.403 


+0 


.4 


-0. 


.7 




15. 


,757 


654 


.757 


-9. 


,4 


+21.0 


.833 


+0 


.5 


+0. 


.6 




16. 


,776 


655 


.776 


+51. 


.5 


-45.9 


.946 


-0 


.5 


+1. 


.1 




17. 


,792 


656 


.792 


+84. 


,4 


-83.3 


13.059 


-0 


.1 


-1. 


,0 




20. 


,846 


659 


.846 


-22. 


.1 


+34.4 


.399 


+0 


.5 


+0. 


.2 




23 


.753 


662 


.753 


-38. 


.3 


+52.7 


.722 


+0 


.3 


+1. 


.2 




25. 


.760 


664 


.760 


+51. 


.7 


-47.7 


.946 


-0 


.1 


-1. 


.0 




26. 


.940 


665 


.940 


+81. 


.4 


-77.6 


14.077 


+0 


.1 


+1. 


,2 




27. 


.836 


666 


.836 


+45. 


.1 


-38.5 


.177 


0. 


.0 


+0. 


.9 


Feb. 


1 


.816 


671 


.816 


-37 


.1 


+50.1 


.731 


+0 


.3 


-0. 


.1 




13. 


.901 


683 


.901 


+82. 


,0 


-79.3 


16.075 


+0 


.3 


0, 


,0 




14. 


,876 


684 


.876 


+43. 


3 


-35.1 


.184 


+1 


.1 


+1. 


.2 




17. 


,838 


687 


.838 


-39 


.3 


+52.1 


.513 


+0 


.1 


-0. 


3 




18 


,834 


688 


.834 


-45. 


.0 


+59.9 


.624 


0, 


.0 


+1. 


.4 




19. 


,798 


689 


.798 


-36 


,9 


+49.4 


.731 


+0 


.4 


-0. 


.7 




20. 


.859 


690 


.859 


-2. 


3 


+13.1 


.849 


+0 


.5 


+0. 


.5 




21. 


.870 


691 


.870 


+60. 


,9 


-58.9 


.962 


-0 


.3 


-2. 


,0 




22. 


.840 


692 


.840 


+82. 


.7 


-80.6 


17.070 


-0 


.1 


-0. 


,1 


Mar. 


14. 


.815 


712 


.815 




+4. 


4 
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